In Search of the Vortex Loop Blowout Transition 
for a type-II Superconductor in a Finite Magnetic Field 



(N 
O 
O 
(N 

o 

Q 



o 
o 

u 

Oh' 



I 

O 

o 



O 

-I— > 



I 

O 

o 



X 
S3 



P. Olsson 

Department of Theoretical Physics, Umea University, 901 87 Umed Sweden 

S. Teitel 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 

(Dated: February 1, 2008) 

The 3D uniformly frustrated XY model is simulated to search for a predicted "vortex loop 
blowout" transition within the vortex line liquid phase of a strongly type-II superconductor in 
an applied magnetic field. Results are shown to strongly depend on the precise scheme used to trace 
out vortex line paths. While we find evidence for a transverse vortex path percolation transition, 
no signal of this transition is found in the specific heat. 
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I. INTRODUCTION 

In pure extreme type-II superconductors, such as the 
high Tc superconductors, the Abrikosov vortex Une lat- 
tice melts via a sharp first order phase transitionn into a 
vortex line liquid as the temperature is increased above a 
critical T™. The properties of this vortex line liquid phase 
have been the subject of considerable investigjation. The- 
oretical argumentsH and early simulationsBolj suggested 
that the vortex line liquid might retain superconducting 
phase coherence parallel to the applied magnetic field, 
within some temperature-interval above Tm- Later, bet- 
ter converged simulationsO found that phase coherence is 
simultaneously lost in all directions upon melting. 

Subsequently, Tesanovi(S has proposed that, for small 
magnetic fields, there still remains a sharp thermody- 
namic phase transition at a temperature r$ within the 
vortex liquid state, associated with diverging fluctuations 
of closed vortex loops, such as drive the superconducting 
transition in the zero magnetic field case. Considering the 
limit of infinite penetration length A, Tesanovic proposed 
that, in a finite field, the fiuctuations of the magnetic field 
induced vortex lines act to screen the interactions of ther- 
mally excited closed vortex loops, in the same way that 
magnetic field fiuctuations screen the vortex loop inter- 
actions of a finite A model in zero applied magnetic field. 
Pursuing this argument, he predicted that the proposed 
vortex loop blowout transition at may be an inverted 
XY transition, as is the case of the zero field Meissner 
transition for the finite A model. 

Following^ Tesanovic's predictions, Sudb0 and co- 
workerajijllj have carried out numerical simulations of 
the three dimensional (3D) uniformly frustrated XY 
model of a type-II superconductor. They claim to find 
evidence for Tesanovic's transition, which they associate 
with the formation of a vortex line path that percolates 
entirely around the system in the direction transverse to 
the magnetic field. 

Most recently, measurementsQ on high purity 
YBa2Cu307 (YBCO) single crystals have found a step 



like anomaly in the specific heat at a temperature high- 
ter than the melting Tm, reminiscent xif an inverted mean 
field transition. It has been arguecLJ that this feature 
may be evidence for Tesanovic's transition r$. 

In order to further investigate this issue, we have car- 
ried out new simulations on the 3D uniformly frustrated 
XY model, both repeating the approach of Sudb0 and 
co-workers, and measuring new quantities that make a 
more direct test of Tesanovic's theory. After correcting 
certain inconsistencies in the earlier numerical work, we 
show that whether or not one finds indications of a vor- 
tex loop blowout transition depends crucially on how one 
chooses to resolve vortex line paths at points where two 
or more lines intersect. Making the choice that favors 
the blowout interpretation, we find the critical exponent 
1/ ~ 1, rather than the value ~ 2/3 expected for an in- 
verted 3D XY transition. Finally, we make high precision 
measurements of the specific heat, in search of a thermo- 
dynamic signature for a blowout transition, but no such 
signature is found. 



II. MODEL 



e is the 3D uniformly frustrated 
^ which models a type-II super- 



The mo 
XY mode^ 

conductor in the limit of infinite magnetic penetration 
length, A oo, and is given by the Hamiltonian, 



i.fi 



cos{0, 



(1) 



Here i are the nodes of a cubic grid of sites, fi = x,y, 
z, are the directions of the grid axes, and the sum is 
over all nearest neighbor bonds of the grid. 9i is the 
phase angle of the superconducting wavefunction on site 
i, Ai^ — (27r/$o) //^^ A • dr is given by the integral of 
the magnetic vector potential A across the bond at site 
i in direction /t, and $o — hc/2e is the flux quantum. 
The argument of the cosine is the gauge invariant phase 
angle difference across the bond. The circulation of the 
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Ai^ around any plaquette of the grid is equal to 27r times 
number of flux quanta of magnetic field penetrating the 
plaquette. We take the magnetic field, B = V x A, 
uniform and parallel to the z axis, with a fixed density of 
flux quanta / — Bo^/^q per plaquette of area a^. The 
couplings we take to model an anisotropic system, 
with Jx = Jy = J±, and Jz ^ J±- 

Simulations are carried out varying the 9i according 
to a usual Monte Carlo scheme; the are held fixed. 
That the do not fiuctuate, and that they give a uni- 
form magnetic field, are the consequences of the X oo 
approximation. Simulations are carried out on x 
cubic grids, using periodic boundary conditions. Except 
where otherwise noted, our runs are typically for 1/4 to 
1/2 X 10^ Monte Carlo passes through the entire lattice. 

While we simulate in the phase angle degrees of free- 
dom, 9i, our interest will be in the behavior of the vor- 
ticity in these phase angles. Let s denote the dual sites 
of the original grid; these are the sites at the centers of 
the unit cells of the grid. Denote by (s,/i) the plaquette 
which is the face of the unit cell centered on dual site s 
with normal in the /t direction, ^ = x,y, z. We define 
the integer vorticity Ug^i piercing plaquette (s, /i) by com- 
puting the circulation of the gauge invariant phase angle 
differences around the plaquette, 

[^'+- - - = 2^(n,^ - fSz^) , (2) 

where the sum is counterclockwise around all bonds form- 
ing the boundary of plaquette (s,/l(), and the gauge in- 
variant phase angle differences are restricted to the inter- 
val [— TT, tt). In a constant magnetic field, the condition 
that the total energy density remains finite can be shown 
to yield the "neutrality" constraint (see section III. A), 

'^nsfj. = fL^LzSzfi , (3) 

s 

i.e. the total vorticity piercing any plane at constant z 
is /L^; these are the magnetic field induced vortex lines. 
The total vorticity in the transverse directions x and y is 
zero. 

Taking the vorticity as the directed bond of the 
dual grid, emanating from site s in direction ft, the vortic- 
ity so defined is divergenceless, forming continuous lines 
that, due to the periodic boundary conditions, must ulti- 
mately close upon themselves. We will label such a closed 
vortex path by the index a, and define the vector Rq as 
the net displacement one travels upon following the path 
a from a given starting point until returning back to that 
point as the line closes back on itself. If Rq ~ 0, then 
the vortex line path is a closed loop of finite extent that 
exists as a thermal fluctuation. If R^a = mLz, m inte- 
ger, then such a vortex line path represents m of the fL\ 
field induced vortex lines; these m lines are mutually con- 
nected to each other via the periodic boundary conditions 
in the z directionB For m > 1, we can say that the m 
field induced lines are geometrically entangled with each 



other. If Rxa = mL^ or Rya = mL±^ then the vortex 
line path winds m times around the system transversely 
to the applied magnetic field. We will be particularly in- 
terested in vortex line paths for which Rza = 0, but R^a 
or Rya ^ 0. Henceforth, we will refer to as the "lines" , 
the set of vortex line paths {a} for which all Rza > 0; 
these are the field induced vortex lines. All other vortex 
paths we will refer to as the "loops" . 

In order to trace out vortex line paths, one needs to 
know how to treat intersections. An intersection is when 
there is more than one vortex line entering and exiting a 
give unit cell of the grid, and it is therefore ambiguous 
which entering segment to connect to which exiting seg- 
ment. It was previously shownES by one of us that the 
method chosen to resolve such intersections can have a 
dramatic effect on the statistics of closed thermally ex- 
cited loops in the zero field f — model. Here, for the 
/ > model, we consider two different schemes, which 
we henceforth refer to as method (i) and method (ii): 

(i) At each intersection we choose randomly, with equal 
probability, which entering segment connects to which 
exiting segment. In the / = model this scheme was. 
found to give results closest to theoretical exjgectations.E^ 

(ii) Motivated by Sudb0 and co-workers ,ffl'E2l we first 
searchllJ through all possible connections to find a path 
a with Rza = and R^a or Rya ^ 0. Such a path winds 
around the system transverse to the field, without ever 
winding around the system parallel to the field. If one 
such path is found, it is selected as a path a' contributing 
to the "loops" , and we then repeat the proceedure applied 
to all remaining vortex paths. When all such transverse 
paths are found, the remaining vortex line intersections 
are resolved randomly, as in method (i). 

Using either method (i) or method (ii) we thus decom- 
pose the vorticity of any given configuration into disjoint 
closed vortex line paths, consisting of a set {a} of "fines" 
and a set {a'} of "loops". 

III. WINDING OF FIELD INDUCED VORTEX 
LINES 

We first attempt a direct test of Tesanovic's theory of 
the T$ transition within the liquid phase. A summary 
of his arguments for the existence of this transition is as 
follows. 



A. Summary of Tesanovic's Theory 

First, a duality transformatioiiEilii0 from the XY 
model of Eq. (0) gives the interaction between vortices 
as, 

nns^] = ^ E [^sf.-fSz^]V^{rs-rs')[ns'^-f6z^] (4) 

where V^(r) is the appropriate anisotropic generalization 
of the Coulomb interaction, with Fourier transform ~ 
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It is this singularity of as (7 ^ that yields the 
constraint of Eq. (|3|) . 

Next, one imagines decomposing the total vorticity of 
the system into lines and loops, 



(5) 



If we define, 



K., EE n^-- - fS,,, (6) 



then J^si^sfj.) = and the Hamiltonian of Eq. can be 
rewritten as. 



n 



s,s' .fj, 



loops 



loops 
s' fl 



(7) 



Tesanovic then argues that a coarse graining of vor- 
tex fluctuations, in the vortex line liquid phase, leads to 
an effective hydrodynamic Hamiltonian on long length 
scales which has the same interaction piece as Eq. , but 
which has a new additive term proportional to '^^^ 6^^. 
The resulting long length scale Hamiltonian then has ex- 
actly the same form as that of a zero field superconduc- 
tor with thermally fluctuating vortex loops, n^°i^P^, and a 
thermally fluctuating magnetic field whose average is 
zero, i.e. the zero field superconductor with a finite pene- 
tration length A. In other words, in this infinite X theory 
at finite magnetic field, the long wave length fluctuations 
of the field induced vortex lines n's""'' screen the interac- 
tion between the vortex loops nl^'JfP*' in exactly the same 
manner as magnetic field fluctuations screen the inter- 
actions between vortex loops in a finite X model at zero 
magnetic field. 

The Meissner transition at Tc in the xero field, finite A, 
model is an inverted 3D XY transition.Ej The high tem- 
perature phase T > Tc has vortex loops on all length 
scales and breaks a globaL|J7(l) symmetry associated 
with a disorder parameter ;ta the low temperature phase 
T < Tc has no vortex loops on sufficiently long length 
scales. The correlation length ^ and rcnopnalizcd mag- 
netic penetration length Aj,; both divergeE3 as ^ \t\~'^ , 
with ~ 2/3 and t EE T - Tc. r-. 

We have earlier carried out numerical simulationst£l of 
this zero field, finite A, Meissner transition. We demon- 
strated that, in this model, magnetic field fiuctuations 
obey the finite size scaling relation. 



F(t,g,L) = {b^iqP)b^i-qiy))/L' 



L-^F{tL^^'',qL,l) 



(8) 



where in the above fi J- i) and h^{qv) = X^s ^ '^'^sp 
is the Fourier transform of of the magnetic fiux density 
00, and q ^ Q, 



hs^i- As L 



F{t, 0,00) 



t<0 
1/C t>0 



(9) 



hence F(t, 0, 00) vanishes below the transition, and in- 
creases continuously from zero as one goes above the 
transition. 



In the present case of a finite magnetic field, if 
Tesanovic's mapping is correct, the Meissner transition 
Tc becomes the transition r$ within the vortex line liq- 
uid phase, and we expect the exact same scaling as in 
Eq. (|^) above, when applied to the quantity bgfj, defined 
in Eq. (^). Taking the hmit of g ^ in Eq. and ap- 
plying to systems with fixed aspect ratio Lz — gL±, we 
expect the scaling. 



(10) 



where t = T - and f{x) = F(a;,0, 1). 
For the directions fi ~ x or y, 



lines 



-'Sfl 



(11) 



is the net vorticity of the magnetic field induced vortex 
lines in the transverse direction fi. The two dimensional 
vector W = {Wx, Wy) defined above is the integer valued 
"winding number" that counts the net number of times 
the field induced vortex lines wind around the system 
in the transverse directions x and y. If {a} is the set 
of vortex line paths that define the field induced vortex 
lines and Ra is the net displacement along path a 

as defined earlier, then X)a -^"J- = WL^, where = 

scaling. 



{Rx,Ry)- We thus expect from Eq. (10) the finite size 



fitL'/n 



(12) 



Note, the neutrality condition of Eq. (H) implies that 
the total transverse vorticity in the system must always 
vanish. For W 7^ 0, it is therefore necessary that any 
such winding of the field induced lines is exactly can- 
celed out by an equal and opposite transverse winding 
of the loops. In the thermodynamic limit, L± 00, 
Eq. (|) implies that (W"^) = for T < Tcs>, and (W^) in- 
creases continuously from zero as one increases T > T^. 
The proposed transition at T$ is thus associated with 
the appearance of infinite transverse loops (see following 
Section IV). 

Another interpretation of the T$ transitkui follows 
from the "two dimensional (2D) boson" modetj of inter- 
acting vortex lines, in which the field induced vortex lines 
are viewed as the world lines of two dimensional bosons 
traveling down the imaginary time axis. For T < T$ 
where {W"^) = 0, the field induced-iiartex lines behave 
like charged two dimensional bosons ,BO with a long range 
retarded Coulomb interaction. In the vortex line liquid, 
T„j < r < r$, where phase coherence is lost parallel to 
the applied magnetic field, the analog 2D bosons are in a 
charged superfiuid state. For T > r$, where {W'^) ^ 0, 
screening by the infinitely large loops n),™'"' results in an 
effective short ranged interaction between the field in- 
duced lines. In this case the winding number squared 
{W'^) is proportional to the superfiuid density of what 
is now an uncharged superfiuid. Thus r$ corresponds 
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to a transition between a charged superfluid and an un- 
charged superfiuid in the analog 2D boson theory. Equiv- 
alently, if one considers the quanta that mediate the in- 
teraction between the analog 2D bosons, the transition is 
from massless quanta for T < T$ to massive quanta for 
T > T$. 

To arrive at Eq. ([l2[) , we considered the transverse 
components of Eq. (|10|). However, we can also consider 
the parallel, fi = z component. Now, 

E = E - JLlL. ^ W.L, , (13) 

S S 

and so we expect the scaling 

(W!) - f.itL't) ■ (14) 

If {a} is the set of vortex line paths that define the field 
induced vortex lines n^g^^^, and Rq is the net displace- 
ment along path a as defined earlier, then J2a — 
{fL\ + Wz)Lz. Thus Wz gives the number of "lines" in 
excess of the average value fL^ set by the applied mag- 
netic field. The neutrality condition of Eq. (^ requires 
that when Wz > 0, there must be an equal and oppo- 
site parallel winding of the loops n!,™'"'. As L — > oo, we 
have = for T < r$, and Wz > ior T > T^. 
Thus a transition at T$ should be characterized by fluc- 
tuations in the number of field induced lines and by the 
appearance of infinite parallel loops directed opposite to 
the direction of the applied magnetic field. 



B. Numerical Results 

To test the above predictions, we have simulated the 
3D uniformly frustrated XY model of Eq. (|^) using a vor- 
tex density / ~ 1/20, anisotropy Jz/J± = 0.02, and as- 
pect ratio Lz/L± = 1, for L± = 10, 20, 30, 40 and 60. For 
these parameters, the vortex lattice melting temperature 
is Tm — 0.24J^, and the zero field critical temperature 
is Tco — 1.14J_L. We compute the transverse winding W 
of the field induced lines, defined by Eq. (pi]), using both 
method (i) and method (ii) to decompose each configu- 
ration into "lines" and "loops" . According to the scaling 
equation (p^, we expect that plots of {W'^) vs. T for dif- 
ferent sizes L± should all intersect at the common point 
t = 0, or T = T^. 

In Fig.|l| we show a semilog plot of {W'^) vs. T/ J± us- 
ing method (i) (random reconnections at intersections) 
for L± = 10, 20, and 30. We see that there is clearly 
no common intersection point of the curves. As Lj_ in- 
creases, {W'^) decreases uniformly over the entire temper- 
ature range. This is in qualitative agreement with earlier 
computations of (W^) by one of us (see Fig. 15 of RefJ^). 
For L± = 60, we have found no net transverse winding 
of the field induced lines at all, i.e. for the length of our 
simulation we had W = 0, for the temperature range 
1.36 < T/J_L < 1.44. 
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FIG. 1: Semilog plot of winding (W^) vs. T/J_i for L_l = 10, 
20, and 30, with vortex density / = 1/20, anisotropy Jz = 
0.02Jx, and aspect ratio = L±. {W'^) is computed using 
method (i), i.e. random connections at intersections. {W'^) 
steadily decreases as L± increases, over the entire tempera- 
ture range. 

Next, in Fig.|^, we show the same quantities but now 
using method (ii) (search first for maximal transverse 
loops), for L± = 10, 20, 30, 40 and 60. We see that as L± 
increases, the curves do seem to approach a common in- 
tersection point, giving aT^ ~ 1.4J^. Note that this T$ 
is above the zero field critical temperature Tco — 1.14J_l. 

From Eq. (p^), we expect that the slopes of these curves 
at r$ should scale with system size as, d{W'^)/dT ^ 

L^l'^ . Fitting each of the curves of {W"^) to a cubic poly- 
nomial in T, we compute their derivatives at the inter- 
section point T = 1.4Jj^, and plot the results vs. L±^ 
in Fig.||. We see that the slopes, to an excellent ap- 
proximation, scale linearly with Lj^, thus suggesting a 
critical exponent v 1. On closer inspection, the data 
in Fig.^ show a small systematic downwards curvature 
about the linear fit; however this curvature can be re- 
moved by assuming a slightly higher critical temperature 
of T = 1.403 J^. Note that this value of ~ 1 is larger 
than the predicted value of 2/3. 

As an alternative method to compute the critical be- 
havior, we can take the scaling equation (p^), expand the 
scaling function f{x) as a polynomial for small x, and do 
a nonlinear fitting to the data to determine the unknown 
polynomial coefficients, T$, and v. To obtain the best 
fit we use a 4th order polynomial and fit only the data 
from the two largest sizes, = 40 and 60. The results 
give T$ ~ 1.403J_L and v ~ 0.96, in agreement with the 
earlier estimates. In Fig. || we show the scaling collapse 
that results from this polynomial fit. There are system- 
atic deviations from the fitted curve on the T > T$ side, 
though these appear to decrease as increases. 

Next we consider the excess parallel winding of the 
field induced lines Wz- As discussed earlier, Tesanovic's 
theory predicts a scahng of (W^) as in Eq. (|l4|). To de- 
termine Wz we count the winding of vortex line paths 
{a} that wind negatively in the z direction, i.e. have a 
net displacement of Raz = —"niaLz, with a positive 
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FIG. 2: Plot of winding {W'^) vs. T/Jx for = 10, 20, 
30, 40, and 60, with vortex density / = 1/20, anisotropy 
Jz = 0.02Jx, and aspect ratio = L^. {W"^) is computed 
using method (ii), i.e. first find all percolating transverse 
loops. Curves of {W'^) intersect at a common point, locating 
T$ ~ \AJ±. Solid lines are guides to the eye only. 
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FIG. 3: Plot of winding slopes d{W^) /dT vs. Lj_ at the 
estimated crossing temperature of Fig.^ T = 1.4Jx, for 
ix = 10, 20, 30, 40, and 60, with vortex density / = 1/20, 
anisotropy Jz = 0.02Jx, and aspect ratio = ix- The solid 
line is the best linear fit to the data. The good fit suggests 
the critical exponent ~ 1. 



integer (R_l may have any value for such paths). Since 
such negative parallel windings must be compensated for 
by excess field "lines", we have Wz = 

However, when we have used either of our tracing 
methods (i) or (ii) , we have never found any such nega- 
tive parallel windings up to the highest temperature we 
have simulated, T = 1.6 Jx. This has motivated us to de- 
fine a third tracing scheme: (iii) we first search through 
all possible connections to find any paths with R^^ < 0. 

In Fig.|^ we show results, using tracing scheme (iii), 
for (Wz) vs. T/J±^ for the same system parameters and 
sizes as used in Fig.|| for {W"^). Note that the values 
of iyV"^) where the curves for different L± intersect are 
exceedingly small. The intersection points appear to de- 
crease in T as L±_ increases, however we do not have the 
accuracy to make any firm conclusions. 

In an attempt to improve the analysis of {W^) we have 
repeated the calculation, but using a new system aspect 




[(T-T^)/J^]Li 



FIG. 4: Scaling collapse of data of Fig. |. {W"^) plotted vs. 
[(T - T^)/Ja_]L^I'", for Lx = 10, 20, 30, 40, and 60, with 
vortex density / = 1/20, anisotropy Jz = 0.02Jx, and aspect 
ratio Lz = L±. Data is fit to a polynomial expansion of 
Eq. (|l|), and ~ 1.4027Jx and u ~ 0.96 determined from 
the fit. Only data from L± = 40 and 60 are used in the fit, 
although data from all sizes are shown in the plot. The solid 
Une is the fitted polynomial curve. 



ratio of = L±/5. This has the effect of increasing the 
value where the curves of (W^) intersect, and so hopefully 
improving our accuracy. We have explicitly checked that 
changing the aspect ratio does not shift the transition 
temperature r$ that is observed in {W'^) (see also Sec- 
tion IV.B). In Fig. I we show results for (W^) vs. T/Jx 
for this new aspect ratio. Again we find no common inter- 
section point for the sizes considered. As L± increases, 
the intersection point continues to decrease. Whether 
this is a failure of the scaling hypothesis of Eq. (|^ , or 
whether we have simply failed to reach the scaling limit 
of sufficiently large {L^ = 24 is the largest value in 
Fig. we cannot be certain. Note that in both Figs.|^ 
and H, (W^) appears to be vanishing at a temperature 
noticeably above the T$ ~ 1.4Jx where the curves of the 
transverse winding, (M^^), intersect. 



We have also tried to fit the data of Fig. g to the scal- 
ing form, {Wz) ~ fz{tLY'^), assuming a non trivial 
anomalous scaling dimension x (although we have no spe- 
cific theoretical reason to propose this form). When we 
do so, we obtain Tc — 1.44, = 0.76, and x — 1.185, how- 
ever our data in the vicinity of this Tc is too scattered 
for us to place much significance on this fit. 

Having used tracing method (iii) to first eliminate all 
possible lines percolating in the negative z direction, we 
can then go and search for all possible transversely perco- 
lating lines and compute the resulting transverse winding 
{W^). When we do this, we find our results for {W"^) vir- 
tually unchanged from tracing method (ii) in the vicinity 
of T$ ~ 1.4Jx. The extremely low number of negative z 
percolating lines at this temperature produces no notice- 
able effect on the transverse tracing. 
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FIG. 5: Plot of z axis winding {W'i) vs. T/ Jx for Lx = 10, 
20, 30, 40, and 60, with vortex density / = 1/20, anisotropy 
Jz ~ 0.02Jx, and aspect ratio Lz = Lx- (W'z) is computed 
using method (iii), i.e. first find all loops that percolate in the 
negative z direction. Solid lines are guides to the eye only. 
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FIG. 6: Plot of z axis winding {Wz) vs. T/Jx for Lx = 20, 30, 
40, 60, 100 and 120, with vortex density / = 1/20, anisotropy 
Jz ~ 0.02Jx, and aspect ratio Lz = I/x/5. {Wz) is computed 
using method (iii), i.e. first find all loops that percolate in the 
negative z direction. Solid lines are guides to the eye only. 



IV. PERCOLATING LOOPS 
A. Summary of Sudb0's Method 

As discussed in the preceding Section III. A, a transi- 
tion at T$ would mark the appearance of infinite trans- 
verse loops, as T is increased. The idea to explicitly look 
for transverse paths that percolate acmss the system was 
first put fogtliJsy Jagla and Balseiro.Eil Later, Sudb0 and 
co-workersBBIlj refined this idea. They defined a quan- 
tity which they denoted Ol, which is the probabihty that 
a vortex path exists which travels completely across the 
system in a direction transverse to the applied magnetic 
field, without ever traveling completely across the sys- 
tem in the direction parallel to the field. If such a path 
exists in a given configuration, that configuration counts 
as unity in the average for Ol\ if not, that configuration 
counts as zero. 

Since having W"^ > in a given configuration neces- 



sarily implies that there is a percolating transverse loop 
in that configuration, there is a close connection between 
the quantities (W^) and Ol- They differ in that (i) for 
a configuration with > 1, and hence with more than 
one percolating transverse loop, the contribution to 
remains unity, rather than increasing with the number 
of percolating transverse loops; and (ii) in a configura- 
tion with two percolating but oppositely oriented trans- 
verse loops, the contribution to Ol will be unity, but 
these loops cancel each other in their contribution to W^, 
which might therefore be zero. 

Since Ol is a pure number one might expect it to be 
a scale invariant quantity, and hence, similar to (W^^), 
plots of Ol vs. T for different system sizes L± should 
have a common intersection point at T$. 

Sudb0 and co-workers' method of searching for such 
percolating transverse paths is siiailar to our method (ii) 
except for one crucial differences They do not require 
that the transverse path closes upon itself; they only re- 
quire that the path start at one end of the system, say 
at a; = 0, and continue until reaching the opposite end, 
X — L±, while keeping the distance traveled along z less 
than Lz, that is the displacement traveled along the path 
a satisfies Rxa = L± and Rza < Lz- Since, by the peri- 
odic boundary conditions, all paths must eventually close 
upon themselves, there are two possibilities for the trans- 
verse percolating paths that Sudb0 and co-workers find. 
We illustrate these in Fig.0: (1) the path a eventually 
closes upon itself without ever traveling the length Lz , in 
which case Rza — 0; or (2) the path a, when followed un- 
til it closes upon itself, does wind up traveling the length 
Lz, with a displacement Rza = Lz] in this case, our 
method (ii) would consider this path as part of the field 
induced vortex lines n'™°'', contributing to the winding 
W, rather than as a transverse loop that contributes to 
■n}°°^''. We wiU call Sudb0's path tracing method (ii') to 
distinguish it from our method (ii). The probability for 
a percolating path using method (ii') we will denote by 
O^; using method (ii) we will denote it by Ol- Paths 
of type (2) will contribute to O^, but not to Ol- We 
will see that there are very dramatic differences between 
these two methods, and that only Ol gives self-consistent 
results. 



B. Numerical Results 

First, we note that if we use method (i) (random con- 
nections) to search for percolating transverse paths, the 
result is essentially the same as found for (W"^) in sec- 
tion III.B. As L± increases, the probability to find a per- 
colating transverse loop steadily decreases for the entire 
temperature range, becoming immeasurably small for our 
biggest system size. Hence we will focus here on methods 
(ii) and (ii'). 

We now consider the computation of O^ using method 
(ii'), the one used by Sudb0 and co-workers, which never 
checks to see how the percolating transverse path closes 
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FIG. 7: Examples of percolating transverse vortex paths. 
Path (1) closes upon itself without ever winding about the 
system in the z direction. Path (2) only closes upon itself 
after winding once in the z direction. Both paths contribute 
to the calculation of by method (ii'). Only path (1) con- 
tributes to the calculation of Ol by method (ii). The thin 
dashed lines represent the periodic boundaries of the system. 



upon itself. We first use parameters Jz = Q.02J± and 
Lz = L±, the same as in section III.B, but a more 
dilute density of vortex lines / = 1/90. For these 
parameters, the vortex lattice melting temperature is 
Tm ~ 0.49J_L, and the zero field critical temperature, as 
before, is Tco ~ 1.14Jj_. These parameters are very close 
to the parameters of one of the cases studied by Sudb0 
and co-workers in Refs. || and ^ (they used / = 1/90, 
Jz = (l/49)Jj_ and ~ L±). Our results for vs 
T/J±, for three system sizes, L± — 30, 60 and 90, are 
shown as the solid symbols on the left hand side of Fig.||. 
These results agree quite closely with those of Sudb0 and 
co-workers (see Fig. 8 of Ref. and seem to show what 
might be a common intersection of the three curves near 
T ~ 0.7J_L. However, we now consider the same parame- 
ters and sizes L± , only using a different system aspect ra- 
tio, Lz = L±/6. The results are shown as the open sym- 
bols on the right hand side of Fig. |[ We sec that there 
no longer appears to be a common intersection point, but 
more importantly, the curves have all shifted dramati- 
cally to higher temperatures. Thus, any value of T$ that 
one might try to extract from depends sensitively on 
the system aspect ratio. We have also considered other 
values of the aspect ratio Lz/L±, not shown here. The 
clear trend is that the sharp rise in shifts to increasing 
temperatures as Lz/L± decreases. But if T$ represents 
a true phase transition, it must be independent of aspect 
ratio. We therefore conclude that and method (ii') 
do not give any self-consistent evidence of the proposed 
vortex loop blowout transition. 

The problems with are even clearer if we consider 
the parameters / = 1/20 and = 0.02Jj^, the same 
ones used for our computation of (W"^) in section IILB. 
In Fig. ^ we show our results for L± = 10, 20, and 40, for 
the two aspect ratios Lz — L± and Lz — L±/2. In both 
cases, there is no common intersection point of the three 
curves for the three sizes, and the curves for the smaller 
aspect ratio are shifted to higher temperatures. Note also 
that, for both aspect ratios, the temperatures at which 
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FIG. 8: Plot of percolation probability vs. T/ Jx for Lx = 
30, 60, and 90, with vortex density / = 1/90, and anisotropy 
Jz = 0.02J_i. Solid symbols on the left are for aspect ratio 
Lz ~ L±; open symbols on the right are for aspect ratio 
Lz — L±/6. O'l^ is computed using method (ii'), following 
Sudb0 and co-workers. Solid lines are guides to the eye only. 
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FIG. 9: Plot of percolation probability vs. T/ Jx for Lx = 
10, 20, and 40, with vortex density / = 1/20, and anisotropy 
Jz = 0.02Jx. Solid symbols on the left are for aspect ratio 
Lz ~ L±; open symbols on the right are for aspect ratio 
Lz — L±/2. O'j^ is computed using method (ii'), following 
Sudb0 and co-workers. Solid lines are guides to the eye only. 



O'j^ rises to unity lie quite significantly below the value 
of T$ ~ 1.4Jj^ found from our analysis of (W^). 

We next consider the the computation of Ol using 
method (ii) (percolating transverse path must close upon 
itself keeping Rza — 0)- In Fig. |l^ we show results us- 
ing the same parameters as were used to compute O^in 
Fig. I i.e. / = 1/90, Jz = 0.02 Jj^, and L± = 30, 60 and 
90, for the same two aspect ratios Lz — Lj_ and Lz — 
L±/6. We see now that for both aspect ratios, curves 
for the three different sizes appear to approach a com- 
mon intersection point, r$ ~ 1.17J_l > Tco ~ 1.14Jj^, 
and that this intersection point is independent of aspect 
ratio (note: for Lz = L±/6, the thinness of the system 
Lz = 5, for L± = 30, presumably makes it too small to 
be in the scaling region, hence it intersects the other two 
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FIG. 10: Plot ol percolation probability Ol vs. T/J± for 
L± = 30, 60, and 90, with vortex density / — 1/90, and 
anisotropy Jz = 0.02J_i. Solid symbols are for aspect ratio 
Z/z = L±; open symbols are for aspect ratio = L±/6. 
Ol is computed using method (ii), where all loops close upon 
themselves. All curves approach a common intersection point, 
Tcj ~ 1.17Jx, independent of the aspect ratio. Solid lines are 
guides to the eye only. 




FIG. 11: Plot of percolation probability Ol vs. T/Ji_ for 
L± = 10, 20, 40, and 60, with vortex density / = 1/20, 
anisotropy Jz — 0.02Jx, and aspect ratio Lz ~ L±. Ol 
is computed using method (ii), where all loops close upon 
themselves. All curves approach a common intersection point 
r$ ~ 1.4Jx, in agreement with the analysis of (W^) in Fig. 
Sohd lines are guides to the eye only. 



curves at somewhat lower temperatures). 

In Fig. ^ we show similar results using the same pa- 
rameters as were used in our computation of {W'^) in 
section III.B, i.e. / = 1/20, J, = 0.02 Jj^ and L,_ = L^. 
We see that the curves of Ol vs. T / J±_ for the different 
system sizes, Lj_ = 10, 20, 30 and 40, all intersect at a 
common point, T$ ~ \AJ^. This is exactly the same 
value as found in our analysis of (VF^) (see Fig.||). 

We therefore conclude that Ol gives a self-consistent 
determination of T$, and that this value is considerably 
larger than estimates one would get from consideration 
of O'^. In fact, estimates of a T$ from all lie helow 
the zero field critical temperature Tco and decrease as / 
increases, while the values determined from Ol all lie 




[(T-T^)]/J^L|' 



FIG. 12: Scaling collapse Ol vs. [(T - T*)/ JxIlY" of data 
from Fig.|l^, for / = 1/90. A reasonably good collapse is 
found for all sizes L±, for both aspect ratios L^/L±, using a 
single value of Tj, ~ 1.168Jx and v = 1. Solid lines are guides 
to the eye only. 



above Tco, and increase as / increases. 

If Ol is indeed a scale invariant quantity, we can pos- 
tulate that it should obey a scaling relation similar to 



I.e., 



(15) 



Based on our analysis of (W'^) in section III.B, we may 
expect v ~ 1. In Fig.|2| we therefore show a scaling 
collapse of the data for / = 1/90 from Fig.|l^, plotting 
Ol vs. [{T — T^)/ J±]L±, where r$ is determined by 
a best fit of the data to the scaling form. We find a 
reasonably good collapse for all sizes, for both aspect 
ratios, using a single value of r$ ~ 1.168J_l. 

In Fig.M we show a similar scaling collapse of the data 
for / = 1/20 from Fig.|l^. Fitting the data for Ol to a 
4th order polynomial expansion of the scaling function, 
we find an excellent collapse, for all system sizes, using 
the parameters T$ ~ 1.399Jx and i' ~ 1.006. These re- 
sults agree very well with the values obtained from the 
scaling analysis of (W^), given in section III.B. The qual- 
ity of the collapse is very much better here than it was 
for {W^) in Fig.|. 

Finally, in analogy with the winding (W^), we have 
also considered the probability Olz to find a vortex path 
percolating through the system in the negative z direc- 
tion, opposite to the applied magnetic field. We expect 
Olz to obey a scaling relation similar to that of Eq. (1^). 
To compute Olz we have used tracing method (iii) in 
which we explicitly search through all possible connec- 
tions to find any such paths. We show our results for 
Olz vs. T/J± for vortex density / = 1/20 and anisotropy 
Jz = 0.02Jj^ in Figs.|l^ and |l^, for system aspect ratios 
Lz — L± and — L±/5 respectively. As with (W^) 
shown in Figs.^ and |^, the intersection points of the 
curves for different sizes appear to decrease in temper- 
ature as L± increases. Again, we cannot say if this is a 
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FIG. 13: Scaling collapse Ol vs. [{T - T.^,) / J a_]L^I'" of data 
from Fig.|l^, for / = 1/20. An excellent collapse is found for 
all sizes L±, using a value of ~ 1.399Jx and u — 1.006. 
These values agree well with those obtained from the scaling 
collapse of {W'^). The solid line is the fitted polynomial curve. 
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FIG. 15: Plot of percolation probability Olz vs. T/J± for 
L± = 20, 30, 40, 60, 80, and 100, with vortex density / = 
1/20, anisotropy Jz ~ 0.02Jx, and aspect ratio Lz = L±/5. 
Olz is computed using tracing method (iii). Solid lines are 
guides to the eye only. 
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FIG. 14: Plot of percolation probability Olz vs. T/J± for 
L± = 10, 20, 30, 40, and 60, with vortex density / = 1/20, 
anisotropy Jz = 0.02Jx, and aspect ratio Lz = L±. Olz is 
computed using tracing method (iii). Solid lines are guides to 
the eye only. 



failure of our scaling hypothesis, or a failure to reach suf- 
ficiently large L^- Also, analogous to our findings for the 
windings {W'^) and (W^), Olz appears to be vanishing 
at a temperature above the r$ = IAJ± where the curves 
of Ol intersect. 



V. SPECIFIC HEAT 

If T$, as determined by (W^) or Ol, does indeed rep- 
resent a true thermodynamic transition, we would expect 
to see some signature of this transition in more conven- 
tional thermodynamic quantities. In the recent experi- 
ments of Ref.pT|, a step like anomaly in the specific heat 
C was observed in the vortex line liquid region, reminis- 
cent of an inverted-*nean field transition. In their numer- 
ical simulations,Q't3 Nguyen and Sudb0 claimed to see 
an anomalous glitch in the specific heat at the tempera- 



ture they identified as r$ from their calculation of O^. 
However, this glitch corresponded to only a single data 
point very slightly displaced above an otherwise smooth 
background; and in the previous section we have demon- 
strated that significantly underestimates T$, hence 
there is no reason to expect any anomaly in C at that 
temperature. 

In this section we report on high precision measure- 
ments of the specific heat C, for the same parameters 
we have studied in the earlier sections. If the T$, as 
found using the vortex path tracing method (ii), is in- 
deed a true thermodynamic phase transition with criti- 
cal exponent ~ 1 (as our scaling analyses found), then 
hyperscaling would suggest a specific heat exponent of 
Q! = 2 — 1. We thus do not expect to see a diverg- 

ing C, however some feature should be present. 

In Fig.|l^ we plot C vs. T/J±, in the near vicinity 
of r$ ~ 1AJ±, for the same parameters / = 1/20, 
Jz = 0.02Jj^, and = L± as used in Figs.|2[ ^ and 
We show results for L_l = 10, 20, 30 and 40, using 
from 3 — 10 X 10^ Monte Carlo passes through the lattice, 
depending on the system size. We find no noticeable fi- 
nite size dependence, and no hint of any feature at all, 
near the previously determined r$ ~ 1AJ±. 

In Fig.|l^ we plot C vs. T/J_l, over a broad tem- 
perature range, for the same parameters / = 1/90 and 
Jz = 0.02Jj_ as used in Figs.|l^ and 12, but for a single 
large system size L± = 30 and = 90. Again we see 
no hint of any anomaly near the previously determined 
T<j, ~ 1.168Ji. 



VI. CONCLUSIONS 

We have carried out detailed Monte Carlo investiga- 
tions of the 3D uniformly frustrated XY model in order 
to search for a proposed "vortex loop blowout" transi- 
tion within the vortex line liquid phase of a pure extreme 



10 



0.70 



0.68 






0.66 








• 


L^= 10 


0.64 










L^ = 20 


0.62 


■ 


= 30 




□ 


L^ = 40 



f= 1/20 
J =0.02 J 

Z 1 

L = L, 



1.36 



1.38 



1.40 1.42 
T/J^ 



1.44 



FIG. 16: Specific heat C vs. T/J± for vortex density / = 
1/20, anisotropy Jz = 0.02Jx, and aspect ratio Lz = Lx, for 
system sizes L±_ = 10, 20, 30 and 40. No hint of any anomaly 
is found near the previously determined T# ~ 1.4Jx. The 
solid line is a guide to the eye only. 
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FIG. 17: Specific heat C vs. T/J± for vortex density / = 
1/90, anisotropy Jz = 0.02Jx, and aspect ratio Lz = 31/x, 
for system size L± — 30. No hint of any anomaly is found 
near the previously determined r<j ~ 1.168Jx- The solid line 
is a guide to the eye only. 



type-II superconductor. Such a transition had been prfe. 
dieted from general theoretical arguments by Tesanovic.Q 
Evidence for such a transition was claimed-|ift, numer- 
ical simulations by Sudb0 and co-workers and in 
specific heat measurements on high purity YBCO single 
crystals.llil We have made explicit measurements of the 
vortex line windings {W'^) and (W^), which are the key 
quantities in Tesanovic's theory. We have re-examined 
Sudb0's calculation of the percolation probability Ol- 

Our results raise several questions concerning 
Tesanovic's theory. We have found that the values 
of {W'^) and (W^) depend sensitively on the precise 
scheme one uses to trace out vortex line paths. For the 
natural choice of random connectivity at vortex line 
intersections, both {W'^) and appear to vanish at 

all temperatures as L —> oo. Only when we specifically 
search first for percolating paths, when computing the 
windings, do we find that the windings converge to 



non zero values above a certain temperature. In this 
case, we find that the transverse winding {W^) obeys 
the finite size scaling form expected from Tesanovic's 
theory, however the critical exponent we find is ~ 1, 
rather than the predicted vxy ^ 2/3 of the inverted 
3D XY transition. For the longitudinal winding (W^) 
we have been unable to find the expected scaling form. 
Whether this is because (W^) does not scale, or because 
our systems are all too small to be in the scaling limit, 
we cannot be certain. It does appear that, upon cooling, 
{W^) vanishes at a temperature above that where (W^) 
vanishes. This would be contrary to Tesanovic's theory. 
However, since we have not succeeded to find scaling for 
(W^), we cannot be certain of knowing exactly where it 
vanishes as L — > oo. 

Independent of Tesanovic's theory, it is natural to 
think that, as temperature and hence vorticity increases, 
the vortex lines may form percolating paths (note how- 
ever that the directedness of the vortex line segments, 
and the condition of divergenceless paths, means that this 
is no ordinary percolation problem!). We have therefore, 
following Sudb0 and co-workers, searched explicitly for 
such percolating paths in the direction transverse to the 
applied magnetic field, as well as in the direction paral- 
lel but opposite to the applied magnetic field. Defining 
transverse percolation as the existence of a vortex line 
path that extends entirely across the system in the di- 
rection transverse to the applied magnetic field without 
simultaneously extending entirely across the system in 
the parallel direction, we have shown that Sudb0's pro- 
cedure, which ignores the transverse periodic boundary 
conditions and does not require the percolating path to 
close upon itself, leads to inconsistent predictions for the 
transition temperature as one varies the system aspect 
ratio. Only by requiring that the transverse percolating 
path close upon itself, without ever winding in the parallel 
direction, do we find a consistent transition temperature 
independent of aspect ratio. The percolation transition 
found this way agrees both in critical temperature T$ 
and exponent with the results from our analysis of the 
transverse winding {W^). We have also computed the 
probability to find a percolating path in the direction 
parallel but opposite to the applied magnetic field. Here, 
analogous to our results for (W^), this negative z perco- 
lation appears to occur at a temperature higher than that 
of the transverse percolation, however we have not suc- 
ceeded to find a clear scaling of this parallel percolation 
probability. 

Note that the transverse percolation transition temper- 
ature T'$(/) that we find increases above the zero field 
transition temperature Tco as the magnetic flux density 
/ increases. This is in striking contrast to the conclusion 
of Sudb0 and co-workers who proposed T$ (/) to decrease 
below Tco as / increases. 

While our results do seem consistent with a well de- 
fined transverse percolation transition, one can ask if this 
is a purely geometrical feature of the vortex line paths, 
or whether it also corresponds to a true thermodynamic 
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phase transition, i.e. something one could detect in a 
suitable thermodynamic derivative of the free energy. To 
investigate this question we have carried out high preci- 
sion Monte Carlo measurements of the specific heat C. 
Our results for C show no feature whatsoever near the 
percolation transition T$, nor do we find any finite size 
effect. In particular we see no evidence for a step like 
feature as was observed experimentally in YBCO. 

To conclude, we have found evidence for a well defined 
transverse percolation temperature within the vortex line 
liquid phase of a model type-II superconductor. The con- 
nection between this transition and Tesanovic's theory 
of a vortex loop "blowout" transition remains unclear. 
It also remains unclear whether or not this percolation 
transition has any observable thermodynamic manifesta- 
tion. 
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